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Abstract 

This study considers implementations of error correction in a sim- 
ulation language on a classical computer. Error correction will be 
necessarily in quantum computing and quantum information. We will 
give some examples of the implementations of some error correction 
codes. These implementations will be made in a more general quantum 
simulation language on a classical computer in the language Mathe- 
matica. The intention of this research is to develop a programming 
language that is able to make simulations of all quantum algorithms 
and error corrections in the same framework. The program code im- 
plemented on a classical computer will provide a connection between 
the mathematical formulation of quantum mechanics and computa- 
tional methods. This gives us a clear uncomplicated language for the 
implementations of algorithms. 

1 Introduction 

The mathematical model of quantum computers is an idealization of a phys- 
ical quantum computer. In a physical quantum computer decoherence with 
the environment causes errors. The use of error correction provides a possibil- 
ity to reduce the effect of errors. A number of mathematical models in error 
correction in the form of an error-correcting code have been developed. We 
will implement some of them in the simulation framework developed by us. 
This framework in Mathematica is a computer language for the simulation of 
quantum computers in classical computers. Within this framework we will 
transform the mathematical model of quantum mechanics into a computa- 
tional code. Thus it will be a straightforward matter to implement quantum 
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algorithms and error correcting codes in this language. More specifically this 
means that it will represent the Dirac notation and theory connected to this 
notation in a natural manner. We build a state in superposition of the com- 
putational basis |0) and |1) act on this state with quantum gates. This state 
will be an increased n-qubit state with the use of the tensor product. Thus 
the n-qubit state will be a superposition of the computational basis |0)®" to 
ll)*^". We will act with one qubit gate or two-qubit controUed-NOT gates 
on this state. Together this will give us a sufficient device for simulating 
quantum computers. 

2 The Simulation Framework 

Let us introduce the part of the program that will be the framework for the 
simulation of error correction. This framework must naturally be part of 
the quantum algorithm that demands protection from error correction. In 
a previous study several well-known quantum algorithms have been imple- 
mented in this framework by the author (see [SI IH E]). We will point out 
that there is a symbolic similarity between our framework and the mathe- 
matical foundation of quantum computing. For this reason we will represent 
the code by a simple modification of Dirac's notation. A quantum state in n 
dimensions can be represented by a linear combination of n numbers of basis 



vectors as {e[0], e[l], . . . e[n]} = {e[0]^", e[0]®"-i ® e[l], . . . , e[n]}. In the 



two-dimensional case a quantum state \(f)) is represented as a superposition 
of two basis vectors, say |0) and |1), known as computational basis (compu- 
tational basis, see [HE])- In this basis a quantum state |0) is represented 



, where a and (5 are complex numbers such as + = 1. We will intro- 
duce some new symbols for the states of the computational basis as follows: 
e[0] = |0) and e[l] = This is the foundation for the structure of the 
program code. For more than one qubit we will use the computational basis 
states e[xi, . . . , a;„] = |a;i . . . where Xj G {0, 1} or the more compact no- 
tation e[y\ = \y), where y = x„2° + ■ ■ ■ + Xi2"~^. We will write the state (p as 
e[0] = ae[0] -|-/9e[l], in analogy to ([H). The operator A acts on the state and 
is usually written as in the quantum mechanical literature. To match 
these symbols we will use the computational symbols A|e[0] for this opera- 
tion. One might regard |xi ... ...?/„)... x„) as |xi ... ?/i ...?/„.. . x„) in 

order to simplify the program code. This will be a computational problem, 
since Mathematica will distinguish between e[xi, . . . , e[yi, . . . , ym], • • • , Xn] 



as 





2 



and e[xi, . . . ,yi, . . . , i/m, • • • , Xn]- The computer must regard these expres- 
sions as equal even if the notations are not identical with each other. As 
an example the expression e[0,e[l], 1] must be equal to e[0, 1, 1] in the code. 
We can bring in the command e[0,e[l],l] := e[0, 1, 1] or the more general 
e[a__, e[6__], c__] := e[a,b,c] to solve this problem. Moreover, the program 
code must be able to handle the linearity of the tensor product. Let e[. ] be 
vectors and a a complex number. We define the tensor product as 



. Two short commands in the program code will implement this definition 
of the tensor product. The command 

e [a , a_. e[x__], b ] := a e[a, x, b] 

will transform e [a] ® ae [x] ® e [c] into ae[a®x^b] = ae[a, x, b]. This command 
is the computational dual to the tensor expression in Dirac's notation \a) ® 
a\x) ® \b) = a\axb). The other command 

e[a___, (a_. e[x__] e[y__]), b___]:= 

^ae[a, x, b]+ ^/3e[a, y, b] 

will transform e[a] ® ^(ae[a;] + P e[y]) ® e[b] to ^ae[a, x, b] + ^l3e[a, y, b]. Let 
U be an arbitrary unitary one-qubit quantum gate. Then U will trans- 
form a one-qubit state e[0], which is represented in the computational basis 
states as e[0] = ae[0] -|- 6e[l], into the state f/| e[0] a(ci e[0] -|- C2 e[l]) -|- 
6(c3e[0] -|- C4e[l]), where a,b,Ci are complex numbers. We add the Mathe- 
matica gate U to the program code as follows: U\ e[0] — ^ ci e[0] + C2e[l] and 
U\ e[l] — s> C3e[0] + C4e[l]. For example, the Hadamard gate H will be added 
in Mathematica as the command iir:={e[0] — > l/A/2(e[0] + e[l]), e[l] — 
1/V2(e[0] - e[l])}. We will define a one-qubit gate as an operator which 
acts on the qubit in position i and leaves the other qubits unchanged. The 
program code must be able to operate with a gate on an arbitrary qubit. 
Consequently, we will define an operator Oi in the Mathematica code. De- 
fined the operator 0^ as Oi = I®'^~^®U (g)J®"~*, which acts on ri-qubits where 
/ is the one-qubit unit operator and U is an arbitrary one-qubit operator. 
Then operator Oi is a function of Oj|e[t'] e[^/']. Similarly, we will define 
Oij as an operator which operates as the two-qubit operator on the qubits in 
positions i, j and leaves the other qubits unchanged. Now we have the tools 
to build the quantum circuit for quantum algorithms and error correction. 



a(e[f] ® e[w]) = («e[t>]) ® e[w\ = e[v] ® (ae[w]) 
(e[fi] + e[t;2]) ® e[w] = e[vi] ® e[w] + e[v2] ® e[w] 
e[v] ® {e[wi\ + e[w2\) = e[v] ® e[wi] + e[v] e[-u;2]. 



(2) 
(3) 
(4) 
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2.1 Quantum Error Correcting 



Quantum computers need some connection to make them controllable and 
will therefore never be completely isolated from their surroundings. The 
quantum computer's surroundings will influence the quantum computer and 
cause errors. The effects of errors caused by inaccuracy and decoherence 
need to be reduced to a minimum. Error correction will help to reduce these 
effects. Define the 3-qubit "logical qubits" denoted (logical-one) and 
|0l) (logical-zero), as it is commonly defined in the literature 

|0i) = |000), |lz.) = |lll). (5) 

Consider some error that will cause a flip on one of the qubits in the logical 
qubits (e.g. |000) change to |010)). This qubit flip will be represented by the 
Pauli operator Xi, where i denotes the position of the flipped qubit. Thus 
we will define Xi = X(S)I®I,X2 = I(g)X(g)I and X2 = I ® I ® X. In 
the same manner a phase flip |1) 1-^ — |1) and the combination of phase flip 
and qubit flip will be represented by the Paul operator Zi and XiZi = —iYi, 
respectively. Moreover, ignore the global phasqj and denote a combination 
of a qubit and phase flip as Fj. A linear combination of the Pauli matrices 
and the identity matrix will then represent an error operator 

Ei = cJi + C2Xi + c^Zi + C4Yi. (6) 

The use of majority voting and a preparation of the state in logical qubits 
([5]) will protect the state against errors that flip a single qubit. This error 
code will not protect against a phase flip. The Shor code [6] will overcome 
this problem. Therefore prepare the state in the logical qubits 

|0,)^(l»°°) + l"'))". |u).(WHm)y\ 

The encoding of the logical states in Shor's code is presented in the left part 
of the circuit in figure [H Some arbitrary error will be simulated by the error 
operator [6] in the middle of the circuit. This operator simulates a one-qubit 
error which can be a phase or a flip error or some combination of these two. 
The decoding is the inverse to encoding, as will be seen in the circuit. 



3 The simulation 

The representations of error correcting code (Shor's code) simulation in 
Mathematica will follow in this section. First deflne the quantum computer 

^Two states which are equal up to the global phase will be equal in the observer's eyes 
(see: p.93]) 
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Figure 1: Circuit for encoding and decoding Shor code 



properties by the code in the hsting 

Listing 1: Definition of register and quantum gates in Mathematica 

e[a ,a_.e[x__] ,b ] :=ae[a,x,b] 

e [a . (a_ . e [x__] . e [y__] ) , b ] : = 

^ae [a , X , b] +^/9e [a , y , b] 

e [a , e [] , b ] : =e [a , b] 

Oj |v_ : =Chop [Expand [v/ . 

(e[x__]:^ReplacePart [e [x] , e [{x} [[i]] ] / . , i] ) ] ] 
Oj ,j |v_ : = Chop [Expand [v/ . [i , j ] ] 
Oj ,j ,A: l'^- ■ = Chop [Expand [v/ . [i , j , k] ] ] 
CN[i_, j_] :={e[x__]:^e[e[x] [[1; ;i-l]] , 

e [SequenceQOMod [e [x] [[i]]+e[[j]] ,2]] , 
e[x] [[i + 1; ; -1]]]} 
T[i,j,k] := {e[x__]:^e[e[x] [[1; ;k-l]] , 

e [SequenceQOMod [e [x] [[k]] + (e[x] [ [ i ] ] * e [ [ j ] ] ) ,2]] , 
e[x] [[k + 1; ; -1]]]} 
H := {e[0]:^l/V2(e[0]+e[l]) ,e[l]:^l/v^(e[0]-e[l])} 
X := {e [0]:^e [1] ,e [l]:^e [0]} 
Y := {e[0]:^-i e[0],e[l]:^i e[l]} 
Z := {e [0] :^e [0] , e [1] :-^-e [1] } 
Id := {} 

This part of the code has defined the register and the quantum gates, and the 
simulation is ready for the quantum circuit. The next part, which expresses 
the quantum circuit, is divided into encoding, simulating errors, decoding 
and measuring. The Enlarge function in the listing 2 will be a simulation of 
the extension of an arbitrary computational one-qubit state to a nine-qubit 
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state. The encoding will affect the Shor code, which is easy to compare with 
the encoding in a quantum circuit [TJ Hence read the Mathematica code from 
inward out. The first quantum gate to implement is the CN4^i gate, the next 
one the CN7 1 gate, and so forth. The simulation of noise is implemented by 
the Error function. After applying some noise we can measure the state and 
return to the initial state. 

Listing 2: Encode and decode 
Enlarge ['?/'_, i_] :=■?/'/. e[x ]:— >-e[x,Sequence@@Table[0,{i}]] 
Encoding [^/'_] : =CN9,7 I (CNgj I (CN6,4 I (CN5,4 I (CNs,! I (CN2,i I (H7 I (H4 I (Hi I 

(CN7,il (CN4,il^))))))))))) 
Error , i] —aiddil^) +bi(Xi\'^) +di(Zi\tlj) +d^(Yi\^) 
Decoding [^/'_] : =17,4,1 I (CN7,i I (CN4,i I I (H7 I (H4 I (Hi I (Tg^gj I 

(CNgj I (CNgj I (T6,5,4 I (CN6,4 I (CN5,4 | (13,2,1 I (CN3,i | (CN2,i 1^) 

))))))))))))) 
Measure [■?/'_] : = FullSimplify [ip / . e [y _ , x__] [y]] 

Let us encode the arbitrary computational initial state tpo and encode by 
the Shor code. The simulation oi a. i : th qubit error on the ^/'i state is 
implemented by the Enlarge function (specific to this example, an error 
in qubit 8). Finally, decode and measure the state. 

Listing 3: Algorithm 

^0 = ae [0]+/3e [1] ; 

ipi = Enlarge [^^o , 9] ; 

ip2 = Encoding lipi^ ; 

ips = Error [ip2 , 8] ; 

^p4 = Decoding lip^} ; 

t/'s = Measure [tp^] 

The output will be (aefO] +/3e[l]) (as + bs + cg — ids), where the global phase 
+ bs + Cs — ids] = lif the error operator is a unitary operator. In fact, the 
global phase can be ignored and the state (ae[0] + /3e[l]) (og + 6g + Cg — idg) 
and (ae[0] + /9e[l]) (ag + 6g + Cg — idg) will be considered as equal in an ob- 
servational point. 
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